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ABSTRACT 



A target must choose a path between some origin and 
destination. The total travel time and the target speed are 
specified, and the target wishes to maximize the "randomness" 
of its track subject to the spatial and temporal constraints. 
Measures of effectiveness are developed against which the 
■randomness" of any path-producing method can reasonably be 
judged. Previous investigations into the scenario are 
reviewed and two models are developed, one using a random 
tour with drift and the other derived from Brownian motion. 
Statistics generated by Monte Carlo simulations for both 
models are compared. While the Brownian motion derived 
process is not always under perfect control of the 
constraints, if the timed arrival constraint may be slightly 
violated then that process performs better against the 
measures of effectiveness and is easier for a target to 
execute than is the random tour with drift. 
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I 



INTRODUCTION 



This thesis will investigate two methods -for constructing 
random target tracks between two speci-fied endpoints. The 
target is constrained to begin its journey at one o-f the 
points and end it at the other a-fter a stated elapsed time. 
Furthermore, either the mean target speed or the maximum and 
minimum speeds will be speci-fied. 

Consider the extreme case in which the distance between 
the two endpoints, the target maximum speed, and the 
speci-fied elapsed time are such that the target is 
constrained to travel directly to its destination. Here 
there can oe r>o randomness to the path, and an attacker needs 
very little target position in-formation in order to make 
per-fect predictions about -future target position. Now 
imagine a less constrained situation in which there is 
considerably excess time. The target is able to travel 
around randoml y during most o-f the period. It might choose 
to move quite randomly at the beginning, suddenly realizing 
at some critical point that it has just enough time remaining 
to travel directly to its destination. However, it is more 
intuitively sensible in this controlled time o-f arrival 
scenario to 'spread* the randomness and e-f-fects o-f the 
constraints evenly across the time period. 
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The two scenarios above give insight about what 
characteristics of target travel are important. In the 
completely constrained case, the target's position does not 
vary at all from its expected position, resulting in an 
easily inferred position throughout the scenario. 
Furthermore, the target always points direcly toward its 
destination, as well as directly away from its origin. For a 
target that wishes to keep its origin and destination secret 
for some tactical reason, completely constrained travel is a 
giveaway. On the other hand, a target that travels about in 
an entirely random manner may never get to its destination. 

It is clear that constrained target travel between two 
endpoints involves many tradeoffs among various constraints 
and choices. This thesis will 1 develop a notion of 
■randomness" by devising measures of effectiveness that 
logically follow from the scenario and that also have strong 
intuitive appeal. Some previous aporoaches to the problem 
will be discussed briefly and then two new 'recipes" for 
target travel will be developed in detail and evaluated 
against the measures of effectiveness adopted. 
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1 1 • MEASURES OF EFFECTIVENESS 



The previous scenario describing totally constrained 
target travel made apparent two dangers to the evasive 
target. First, the target risks divulging both its origin 
and destination because it always points toward its 
destination and away -from its origin. It -follows that the 
less constrained a target is, the less it necessarily points 
to destination and away from origin. Exactly how much course 
freedom is gained as the target becomes less constrained wi 1 1 
be investigated in detail later. One might counter the 
importance of this pointing by claiming that, at best, the 
attacker obtains only a line of bearing to the origin or 
destination. While true, this argument neglects the 
possibility that other targets may pass the same way, with 
the same origin or destination. Sooner or later the attacker 
will get lines of bearing that cross with a regularity 
sufficient to specify the critical positions. Since origins 
and destinations might be important enough to be kept secret, 
one reasonable measure of effectiveness against which to 
evaluate any set of target paths is the absolute angle 
between the target's present course and the course from the 
target's position to its destination sampled at specified 
time intervals <c.f., Figure la). For the totally 
constrained target this angle will always be zero. In order 
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to compare two different path generation procedures it will 
be help-ful to produce a plot o-f this measure against time if, 
in fact, it is time dependent. Otherwise, it will be useful 
to calculate the mean pointing angle, its standard deviation, 
and perhaps an empirical density function. 

The second undesirable quality of a completely 

constrained path is that the target's position always falls 
exactly on its baseline position, where baseline position 
is defined <is the point where the target would be at any 
time of its Journey if it travelled directly to its 

destination at a constant speed. For example, if the distance 
between origin nnd destination is sixty distance units and 
the specified elapsed time for the problem is thirty time 
units, then the rate of travel along the baseline is two 
distance units per time unit. Ten time units after the start 
of the journey, the target's position is on the baseline 
twenty distance units from the origin and forty distance 
units from the destination. The baseline is important 

because the more constrained the target is, the closer it 
must stay on the average to the baseline. Increasing 
constraint on the target by decreasing allotted travel time 
or by increasing baseline length straightens out the target 
path with the result that inferred positions are easier to 
obtain. Very importantly, even if the target is not forced 
to stay near the baseline as a result of being constrained, 
but chooses to do so by hovering or zig-zagging along the 
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line, an inferred course and position are easily obtained. 

Even worse, such an interred course will most likely point to 

the destination or from the origin. Accordingly, the second 

measure of effectiveness shall be the mean square radial 

distance between the present target position and the baseline 

target position at specified time intervals during a journey 

2 

(c.f., Figure lb). This statistic will be denoted as ECR^. ]. 
The squared radial distance has been chosen rather than the 
radial distance because previous investigations into the 
problem, which will be cited in the next chapter, have tended 
to use radial distance squared. It is easier to handle 

analytically than is R^, and can be viewed as the sum of the 
x and y component squared distances from the component 
baseline positions at any time t: ♦ 



R 2 = <X.-£ J 2 -»• .) 2 < 1) 

t t xt t yt 

2 

It would be not difficult to in':egrsite R^ over the whole 

time interval for an individual path in order to get a scalar 

value for this performance measure. Likewise, the mean 

2 

square radial distance ( ECR.. 1 ) could be integrated for 

each particular path-producing recipe in order to obtain a 

scalar measure for that method. Then, two different path 

methods might be compared and, all other measures being 

equal, the method with the highest value selected. However, 
the reduction of this measure of effectiveness to one scalar 
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can result in great loss of information. Two path generation 

procedures might have the same value by the scalar measure, 

but be charac ter i st i cal 1 y different. One might tend to have 

relatively large values of EtR^. 3 near the origin and 

destination, while the other method might have small values 

there, with large values for t in the middle of the travel 

period. Furthermore, each of these two different path 

recipes might be more desirable than the other under 

different circumstances. For instance, if it is not 

important to keep secret the origin and destination because 

they are already known to the enemy, but it is important to 

be as undetectable as possible <or if detected, as difficult 

as possible to redect) in between origin and destination, 

then the target would probably prefer the latter of the tvo 

path-producing methods. In any case, for any path-pr oduc i ng 

recipe, the captain of the target should be able to look at 

2 

the plot of EIR^ 3 as function of time. Thus, the measure of 

2 

effectiveness itself shall be either a plot of EtR^l against 

time, the function that describes i t as such, or a list of 

ordered pairs with time as the first element and mean square 

radial distance as the second. 

A careful observer might object to this measure of 

effectiveness. While it is true that a highly constrained 

2 

target will exhibit a small EIR^ 3 for all t and have easily 

inferred positions, this does not necessarily mean that a 
2 

large EtR^ 3 will guarantee that target position will be 
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diff icul t to infer 



2 

In fact, if maximizing EIR t 1 over the 

entire time period is taken as the only measure of 

©■f "f ec t i veness , it should be almost as easy to infer future 

position from a few past positions as it is in the totally 

constrained case. Imagine a baseline ten distance units 

long, a maximum target speed of four distance units per time 

unit (4 d/t) and a specified travel time of four time units. 

It seems sensible that the target would travel at maximum 

speed the entire time in order stay as far away from its 

expected position as possible. For a good analogy, imagine a 

stick ten distance units long with a piece of string sixteen 

distance units long attached to it, each end to each end. 

Then, the target might choose a track to maximze the area 

between the stick and the string. <The analogy is not 

perfect however; in the stick and the string example the 

radial distance <R^> integrated over the entire . time period 

is being maximized rather than the radial distance squared 
2 

(Rj. ). One possible path choice is to drive the target 

straight away from the baseline at about a fifty degree angle 

for eight distance units, then turn back and go directly to 

destination. The area of the isosceles triangle thus formed 

with the stick as the base is approximately 31 distance units 

squared. But a clever hard working target could drive in a 

circular arc, thus enclosing approx imatel y 39 distance units 

2 

squared. This path, while producing a large EIR^ I at each 
time t, is not very "random" and might be easily targeted. 
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The third and final measure of effectiveness is important 

in at least two ways. The first pertains directly to the 

shortcoming of the second measure of effectiveness taken by 

itself. Given a certain recipe for producing paths, it is 

necessary to produce many of them in order to obtain 

statistical estimations of both of the two measures of 

effectiveness already chosen. It may be possible to apply 

seme non-linear regression method to any path, and to extract 

an estimate of future position to use as the expected 

position for any time t, rather than use baseline position as 

already defined. Such a method would quickly identify the 

smooth curves produced by any path recipe which maximized the 

mean square radial distance only. But the task would be quite 

difficult and costly, especially when applied to the several 

hundred paths necessary to produce good statistics. In order 

to save effort and money, the third measure of effectiveness 

will subjectively judge how represen tat i ve paths from each 

generating method "look*. Against this measure, any path 

2 

which maximizes ECR^ ] alone will be rejected at a glance if 

it exhibits long straight legs or a predictably curved path. 

However, if a path does not look so regular as to be 

2 

predictable, then a large ECR t ] is desirable because it 
means the target is staying away from the critical line 
between origin and destination. 

A second very important reason that paths will be 
visually inspected is that they must be able to be executed 
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by an actual target. Hence, an experienced person must Judge 

whether or not a path is practical regardless of its 

statistics. For example, a high -curvature path which 
2 

maximizes ECR^ 3 may not be impl emen tabl e. 

There are, then, three measures of effectiveness that 

shall be used to judge path-producing recipes. Though they 

were introduced in a different order to facilitate logical 

development of them, they will be applied f.o a path procedure 

as follows. First, a representat ive path must "look good" by 

being practical, executable, and random looking (no 

noticeable regul ari ties) . If a path can pass this first 

important test, then several hundred will be generated using 

the same procedure, and the mean angle (and its standard 

deviation) between present course and cojrse from present 

position to destination will be calculated at specified time 

intervals. Then, the mean square radial distance between 

present position and baseline position will be calculated at 

2 

each specified time interval. The ELR^ 3 will thus be 

estimated at various stages of path completion. If two path 

generating recipes both produce paths that meet the origin, 

destination, time, and speed constraints, and "look good" by 

the first measure of effectiveness, then the path which 

exhibits the least pointing to destination as judged by the 

2 

second measure of effectiveness and the greatest ECR^ 3 at a 
given stage of path completion shall be judged to be the more 
desirable path. Some comparisons will no doubt result in a 
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situation in which both procedures produce paths which "look 
good", but each of the two methods has a better evaluation 
than the other in one o-f the other two measures. One might 
try to -form some weighted combination o-f the two measures, 
but this is dangerous because measures o-f e-f -f ec t iveness do 
not combine well; at minimum, their units are not generally 
on the same interval scales. Such reductionism is not 
necessary anyway. One need only regard both o-f the last two 
measures and decide when it is advantageous to weight one 
subjectively over the other. In fact, any formal weighting 
system would probably not be able to capture all scenario 
dependencies as well as a subjective weighting. 
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III. PREVIOUS INVESTIGATIONS 



While several approaches to the problem of randomizing 
target motion are possible, one o-f the most appealing is 
investigated in detail by Washburn [Ref. 13. In his model a 
target takes a random tour by choosing its direction o-f 
travel from a uniform <0,2iO probability distribution and its 
length of travel on the selected heading from an exponential 
distribution with parameter X (mean number of turns per time 
unit). Thus, turning points are the jumps of a Poisson 
process with parameter X, and at any point during the process 
both the backward and forward recurrence times are themselves 
distributed exponentially with parameter X. This property of 
'memory 1 essness' is very appealing. The probability that the 
target does not turn by time < t+a) , given that it has not 
turned by time t is the same as the probability that the 
target will not turn during the time interval <9, a). In 
other words, in deciding at any point in time when to turn in 
the future, the target does not remember how long it has 
already traveled on the same course. Hence, an attacker may 
not infer either that the target will stay on a new course if 
it just turned to it, or that the target will turn soon 
because it has been on the same course for a long time. 
Because of the desirability of giving an attacker so little 
information, each leg of target travel is chosen by draws 
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■from the -fol lowing distributions: 



•f (O) = 



f<t) = 




f 

< 

V 



Xe 



-xt 



0 



0 < e < 2« 

otherwise; 



t > 0 

otherwi se. 



< 2 ) 



After much math«?mat i cal manipulation, Washburn derives 
the probabi 1 ty density o-f the target's radial distance -from 
its origin, given no initial course information, to be: 



-f <r , t)=t l/2*<Vt> 2 JCXt/< l-r 2 )^expC-Xt< l-< l-r 2 >*> ] . (3) 



He also makes some interesting observations: 



1. The larger X is the more the distribution piles up 
around the origin. 

2. Given no turns, the target is uniformly distributed on 
a circle of radius Vt (throughout this thesis, a 'V' 
will denote a scalar which is the magnitude of the 
vector denoted by 'V'). 

3. Very remarkably, given that two steps have been 
completed at time t, the target is uniformly 
distributed in the circle of radius Vt , not including 
the c i rcumf erence. 

The first two observations go hand in hand. Recall from 
Chapter I that a target which hovers about the baseline or 
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zig-zags along it has a small E[R + ~] since it is rarely -far 

■from its baseline position. Washburn's distribution 

quanti-fies exactly how much the target remains near its 

expected position relative to how much the target "hovers", 

as indicated by the parameter lambda. The only di-f-ference is 

that Washburn's model has no dri-ft, resulting in the expected 

position being a single point rather than a baseline o-f many 

points as in the drift case. A turning parameter of lambda 

equal to infinity creates the degenerate case in which the 

target is distributed on top of its origin with probability 

equal to one. The other degenerate case occurs in 

observation two when the target does not turn and is then 

necessarily on the circle of radius Vt. Just as noted in 

2 

Chapter I, E[R^ ] may be high in this form of degeneracy, but 
given any two bits of position information, future position 
can be inferred perfectly because the target has been 
travelling in a straight line. 

Though it is interesting that the target is uniformly 
distributed on the disk of radius v*t if it just happens to 
have finished leg two by time t, this fact really does not 
help a target evade an attacker. Given the set of 
instructions by which paths are generated, it would be purely 
coincidental that the second leg is completed right when the 
attacker looked for the target < the probabi 1 i ty is zero). 
One might say that the target should plan to finish the 
second leg just when the attacker is expected. For example, 
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the target might know when it conducts some type of evolution 
that makes itself more detectable. Also, it also might know 
from how far away the attacker must come and at what speed, 
so it has a good idea when its window of vulnerability is. 
But if it uses this information to plan two legs then it is 
not following a random tour as prescribed by the probability 
distributions set forth. If the target is going to break the 
rules, it may as well just pick a point at random inside the 
circle of radius Vt and head for it. This procedure will 
guarantee that the target has sampled uniformly over the 
di sk . 

Belkin [Ref. 23 does pioneering work in comparing 
Washburn's random tour process with a Gauss-Markov process 
for describing diffusion, the Ornstei n-Uh 1 enbeck < I QU> 
process. As a result, Belkin finds that the mean square 
radial distance of a target from its expected position at 
time t in a random tour process is: 





ECR t 3 = 2<V /X > < e 


+ Xt - 1) 


• 


( 4) 


Belkin [Ref. 


33 further embellishes 


his analytic work 


on the 


random tour 


process by deriving 


the mean 


square 


radi al 


distance for 


a random tour wi th 


arbi trary 


course 


change 



distribution to be: 

ECR. 2 ] = 2<V 2 /X 2 > < 1 - |A 2 - 1* 2 > <e“ Xt + Xt - 1), < 5) 

t w x y 
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where , 




cos© dF<©> 



e 



( 6 ) 



2 * 



y 



/ 

0 



sin© dF<©) 



Here F<©) is the cumulative distribution -function o-f © -for 
which there may or may or not be a proper density function. 
The subscript w on Y indicates that the velocity is total 



velocity through the water, comprised of both a drift 
component and a randomizing component which will be explained 
later. . Notice that for the uniform distribution of the 



results as expected. 

The notion of an arbi trary course change distribution is 
important because the expected position of the pure random 
tour process is the origin at time zero. As such, the 
process will never cause a target to migrate toward its 
destination. Therefore, it is necessary to introduce some 
bias into the selection of the courses in order to keep the 
target moving in the correct direction toward its 
destination; and even a course distribution alone will not 
cause a target to visit the destination without the addition 
of some further constraining process. While an infinite 



Washburn's random tour model, H = F =9, and equation <3> 

j x y ’ ^ 
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number of arbitrary course selection distributions are 
possible, several look promising when applied to the present 
problem. Loan© CRef. 4] suggests: 

f<6) = i/2 X (l + occose) for 0 « <0,210 <7> 

a « <9, 1) , 

where bias in the positive x direction can be controlled by 
the selection cf the parameter a. Because constraining 
target travel by course distribution alone only causes the 
target to migrate toward its destination, without giving any 
guarantee that the endpoint will be visited, some process 
must be devised to meet the visitation constraint. One of 
many possible course distributions, which also injects the 
time constraint into course selection, is: 

0 normal ECUS, Arcos<T ^/T^) ] (8) 

where C'JS is the course to the destination at the beginning 

of the leg, is the time it would take to go directly to 

the destination at mean speed or maximum speed (whichever is 

selected as important for a particular run), and is the 

time remaining in the problem. Notice that in the most 

constrained case in which T , = T , the variance of the course 

a r 

distribution equals zero and the target travels straight to 
its destination. While this course distribution causes the 
target to travel over the destination, it does not control 
when it will cross it. Hence, it does not solve the 
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controlled time of arrival scenario devised in Chapter I. 
Belkin [Ref. 53 meets the visitation constraint with the IOU 
process and is able to distribute the contraints evenly 
across the whole path through optimal control of velocity. 
This means that the target arrives at destination when it has 
to, without either traveling to i t strai gh taway and hovering 
for the remainder of the period, or by randomizing travel 
until the last possible opportunity to get straight to the 
destination on time. Nonetheless, mean square radial 
distance from baseline position will necessarily need to 
begin decreasing at some point in order to guarantee that the 
target visits its destination. Two observations made by 
Belkin motivate the approaches taken in this thesis for 
devising path recipes. He states ERef. 6 3 that "...it is 
possible to approximate a random tour process with arbitrary 
course change distribution by an IOU process with linear 
drift." He also notes [Ref. 73 that “...as the constrained 
process approaches the terminal constraint at T, the process 
behaves precisely like an unconstrained process running 
backwards in time from T." T, in this instance, is the total 
time alloted for the problem and the terminal constraint is 
the requirement for the target to be at the destination at T. 

As a result of Belkin's insights, the first of two path 
producing recipes that will be developed in this thesis will 
treat constrained travel between two endpoints as a random 
tour process with drift (distance between endpoints divided 
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by total tilne) , executed from both endpoints with an attempt 
to connect the two separate paths somewhere in the middle. 
The drift vector applied to the path beginning at the origin 
will "blow" the target toward the destination, while the 
drift for the path starting at the destination will "blow" in 
the opposite direction, toward the origin. This approach 
constitutes a discrete approximation of the IOU process with 
drift, wherein the random tour out of the destination is one 
half of the total problem run in reverse. 

Belkin's solution to the problem of constrained target 

travel is the IOU process and he has developed a computer- 

simulation named IOUTRK for which he presents some sample 

paths CRef. 83. While the genesis leading to his adoption of 

the IOU process begins with Washburn's random tcur , Belkin 

proceeds to attempt an approximation of the random tour 

process using Brownian diffusion [Refs. 9,103. However, his 

investigations lead him [Ref. 113 to the very clever 

realization that the functional form of the mean square 

2 

radial distance, E[R^ 3, is exactly the same for the IOU 
process as for the random tour after making only two simple 
parameter substitutions [Ref. 123. He concludes his analysis 
[Ref. 133 by stating that "...if one is constrained for 
theoretical or computational reasons to approximate the 
motion of a randomly touring target by Gaussian diffusion, 
then the Ornstei n-Uh 1 enbeck displacement model is to be 
preferred to the Brownian motion model ." 
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The random tour process is important and has many 
desirable properties that will be discussed -further when the 
■first new path recipe in this thesis is developed and 
analyzed. As indicated previously, the random tour will not 
be approximated by any other process; it will itsel-f be 
executed in a new way. The second model to be o-f-fered in 
this thesis will examine Brownian motion, not as an 
approximation of a random tour, but as a basis for solving 
the constrained target motion problem in another way. Though 
rejected by Belkin as an approximation of the random tour 
process, Brownian motion stands alone as a method for solving 
the problem at hand in a very novel way. 
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IV. PROCEDURE 



A. RANDOM TOUR WITH DRIFT 

In order to evaluate the random tour with drift against 
the measures of effectiveness delineated in Chapter II, a 
Monte Carlo simulation i i* devised. Specifically, two random 
tour processes are executed, one from each endpoint, in order 
to guarantee that the destination (actually the starting 
point of the second process) is visited as required by the 
constraints. Recall from Chapter III that the reason for 
executing two separate processes stems from Belkin's 
observation that the path approaching the destination looks 
like an unconstrained random tour run in reverse. Hence, the 
path which visits the destination is constructed by executing 
an unconstrained random tour originating at the destination. 
For all paths, and without loss of generality, the origin is 
zero on the x axis of a Cartesian coordinate system in two- 
space and the destination is however many distance units 
desired in the positive x direction. The random tour 
beginning at the origin is “blown" in the positive x 
direction by the drift vector and the process originating at 
the destination is “blown“ in the negative x direction by an 
opposing drift vector of the same magnitude. The goal is to 
generate two paths that “blow* into each other somewhere 
between origin and destination. At the outset, the distance 
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between origin and destination, the time period, the mean 
target speed through the water (V^) , and the target turning 
rate < 1 ambda) must all be specified. Drift velocity is then 

in the positive or negative x direction, depending on which 
of the two random tours is being executed, with a magnitude 
equal to total baseline distance divided by total problem 
time. Figure 2 below illustrates the other important 
velocity, the randomizing velocity V^ , the magnitude of which 
must be computed before the random tour can begin. 

Figure 2: 

Water Velocity as a Composition of Drift 
and Randomizing Velocities 



♦ 




Vft sin © 



Because the target must have a mean speed through the water 
of V , and V. . is determined by baseline length and total 
run time, there is no freedom in choosing mean V^ . Recall 
that the angle theta is the direction of target travel chosen 

from the uniform, probabilty distribution <0,2lO and notice in 

^ 

Figure 2 that V = V , + V . 

w d r 
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Then : 



2 * 



E IK> 2 ] = EC IV +v : 2 ] = 
w dr 



1/2KCM 2 sin 2 0+(U cos^+U J 2 3 
r r d 



2 * 



1/2* I <V 2 + M 2 + 2V M .cose) de 
r d r d 



< 9) 



8 




+ 






2 



d 



Hence, given V . and a desired mean V , the randomizing 

d w 

velocity for the random tour with drift should be: 

U = <U 2 - V 2 ) * < 10) 

r w d 

Notice that as the drift speed approaches the water speed, 
the randomizing speed goes to zero. This is reasonable 
because when the target has just enough time at a given speed 
to go straight to the destination, all its velocity goes 
toward matching the drift velocity and there is no excess 
velocity left over for any randomization; the target is 
totally constrained. 

With these prel imi nar i es taken care of, the simulation is 
executed beginning with the "left" path <out of the origin at 
zero on the x axis) according to the following procedure: 



1. Choose the direction of travel from U<0,2*>. 
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. Select the length of each leg of the process as U t, 
where t is determined by a draw from an exponential 
distribution with parmeter lambda. 

3. Compute the actual vector V) describing travel for the 

w 

leg by adding a prorated drift vector V , the magnitude 
of which equals [<t / total problem time) x baseline 
length], to which was determined in steps 1 and 2. 

4. Continue executing steps 1-3 until the cumulative 
travel time -for all left side legs is at least .475 of 
the total time -for the whole problem. 

After one left side path is simulated, a similar process is 
executed from the “right" side (destination) using an 
opposite drift direction. As a result there are now two 

random tour paths which migrate toward each other at the 
determined drift rate. The two processes use up at least .95 
of the total time alloted for the target travel from left to 
right. Now, a circle with radius .05 x total time x desired 
water speed is drawn around the left endpoint. The 
proportion of time left in which to join the two sides is not 
arbtrary, but depends on the length of the total time of the 
problem. Here, the total time is assumed to be between 30 

and 60 time units, making the joining period between 3 and 6 
time units. During this connecting period the target should 
maneuver as suggested by X, laying down a path that will make 
the connection. If the right path endpoint falls within the 
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circle, then a "match" has occurred and the two connected 
paths constitute a complete sample path satisfying all 
constraints. If a match is not realized in a specified 
number of iterations of the right path, then a new left path 
is executed and another attempt to generate a matching right 
path is made. The number of right paths that are generated in 
attempting to match a given left path and the number of left 
paths that are generated after a right path match iteration 
limit is reached can both be controlled by parameters in the 
simulation. The entire process is continued until 100 
matches < or any other number desired) are obtained or until 
the limits for the number of attempts are reached. The 
motivation for this procedure of matching is twofold. First, 
desired statistics can be generated to test the symmetry of 
both sides since many right side and left side paths will be 
generated during the quest for 108 matches. Secondly, in 
attempting to get matches., the number of right sides allowed 
in order to match a given left side can be strictly 
controlled. The lower the limit on right sides, the more 
likely i t wi 1 1 be necessary to generate another left side. 
In this manner, the total number of left sides generated in 
order to obtain 188 matches can be much greater than 180, and 
the number of right sides will be even greater yet since 
right sides are constrained to match left sides and not vice 
versa. This path generating method results in two complete 
distributions of paths for the left side. Some left hand 
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paths result in matching right hand paths, while other left 
paths are rejected because no right paths can be generated 
(within limits) to match them. The retained left paths 
comprise a set which will be called the "clipped* set, while 
all paths, rejected or not, will be called the "unclipped* 
set. These names are reasonable because the paths that match 
are more constrained; they have not been rejected as 
unacceptable. A comparison of the statistics of interest for 
both left hand sets will indicate how much degradation 
results from constraining the clipped set. Additionally, 
there is a third set of paths which contains all the right 
paths which match left paths. Statistics are generated for 
this set in order to check symmetry with the left hand 
cl ipped set . 

In order to facilitate the generation of valid statistics 
it is necessary to sample the left and right paths at the 
same specific time intervals. While it is much easier to 
sample a path at the end of each leg, the resulting 
statistics do not provide consistent comparisons along the 
time line. Accordingly, statistics are generated at times 
determined as a function of the turning rate, so that the 
probability of generating statistics twice on the same leg 
is .801. The statistics generated at each time check for all 
three sets of paths are: 

1. Mean magnitude of the angle between present course and 
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course to destination -from present position; variance 
and standard deviation o-f that mean angle. 

2. Mean square radial distance between baseline position 
(pure drift prorated for cumulative run time) and 
actual position. 

Using the procedure described here, statistics are generated 
for a water speed of 4.0 dist/time and drift rates of 3.75, 
3.0, 2.0, 1.0 dist/time, all for turn rates of 1.0, 2.0, and 

4.0 per time unit. There are a total of twelve cases. 

B. BRQWNIAN-DERIMED MOTION 

Since Robert Brown first observed the highly irregular 
motion of suspended pollen particles in 1827, mathematicians 
and scientists have spent a great amount of time and effort 
investigating Brownian motion, and its widespread 
applicability to naturally occurring events makes it a 
reasonable candidate for randomizing target motion. From a 
probaba 1 i st i c viewpoint Freedman [Ref. 14] defines normalized 
Brownian motion to be a stochastic process (B(t):0 < t < «> 

on a sample space 0 with properties <a-c) for points (0 f 0: 



(a) 


B(0,O)) =0 for each 0), 




< b) 


B(.,0)> is continuous for 


each 0), 


<c) 


for 8 < t . < t 0 <. . .< t 
12 n 


. < t , the increments B(t.), 
-In’ 1 




B< t 2 > - B< t j) , . . . , B( t R > 


- B(t^_^> are independent and 



35 



normally distributed, with means 0 and variances t^, 

*2 " V***’ t n ” t n-r 

Additionally, this process is Markovian, meaning that the 
■future is conditionally independent o-f the past, given the 
presen t ; 



where the 
var i abl e X 
Once 

shal 1 be 
coordi nate 
component 
coordi nate 



pc x t+1 < i t+1 ! X g - i 0 ,..., x t - i t > 



p< X t + 1 5 A t + 1 ! x t - V 



< 1 1 ) 



i's are elements o-f the state space -for the random 



again, the problem o-f constrained target motion 
-framed in two-space using a standard Cartesian 
system. Taking position to be specified by each 
separately, the conditional probabi 1 i t i es -for each 
at time t, given the initial constraints are: 

PC X<t> = x ! X<0) = 9, XCT) = L> 



< 12) 



PC Y< t> = y i Y< 0) = 0, YCT) = 0} 



where t is cumulative run time, T is total alloted time and 
(0 < t i T) , L is the positive x coordinate o-f the 
destination, and 0 is the y coordinate o-f destination. 
Hence, the Cartesian layout is exactly as it is -for the 
random tour case. Conditioning the Brownian motion 
constructs a “Brownian bridge 0 between the origin and the 
destination, by which the initial problem constraints are 
met. The x component conditional probabi 1 i ty expands to: 
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P < XCt) = X i XO) - 9 > . P { X<T) = L ! X < t ) = x > 

P ( X<T) = L ! X<0> = 0 } 



= P < X<t) - x. X'lT) = L ! X(0) = 0 > 
P C X<T) = L i X<0) = 0 > 



< 13 ) 



and the conditional density when X<t) is distributed Gaussian 
becomes : 



exp<-l/2£ <x-0)/< dc'^S ] 2 > 



exp<-l/2£ <L-x)/<d<T-t) *) 3 2 > 



(210 



<210 



<T-t) 



< 14) 



exp<“ 1/2C < L-0) /< <?t^) ] 2 > 



<210 



<? . T 



exp< - 1/2 Q) . T* 

it x X 

<210 ‘ . t . <T — tp 



where : 

Q = <x 2 /d 2 ) * i < L-x) 2 /C O 2 <T— t) ] > - <L 2 /d 2 T). 

A-fter completing the square -for Q and simpl i-fying, the 
conditional density reduces to: 

exp<- 1/2C <x-Lt/t)/<d 2 t< 1-t/T) ) ] 2 > 

<15) 

<210^ . <S . <t£ l-<t/T) ]> . 

Hence, X< t) and Y< t> £by a similar derivation] are distributed 
X< t) ' Normal £ Lt/T, d 2 t<l - t/T) ] 

< 16 ) 

Y<t)^ Normal £ 0, <? t< 1 - t/T)]. 
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In both distributions <r is a physical parameter that may be 
specified, and the total variance describes a parabola when 
plotted against t. The maximum variance is (</“)< T/4) and 
occurs at t = T/2. Very importantly, the variance is zero at 
both the origin < t = 0 ) and the destination ( t = T ) . 

One of the measures of effectiveness chosen in this 
thesis against which to judge a target path examines how a 
path "looks" and whether it is executable by a crew. 
Brownian motion, as defined, is not executable because it is 
a continuous process with an everchanging velocity, and as 
such, is impossible for a large target to duplicate. However, 
there is no reason why the continuous process cannot be 
sampled at various times and the selected points be made the 
endpoints for legs of straight line travel. Of course 1 , as the 
time interval between samples is lengthened, the* linear 
approximation connecting the sample points becomes less 
"Brownian". Nonetheless, the process still retains vestiges 
of its Gaussian properties, and the flavor* of Brownian 
mot i on . 

An exponential distribution should be used to determine 
the sample times for the same reason that is advanced for the 
random tour: the memoryless property applied to course change 
is advantageous to the target. Exactly how to sample from 
the normal distributions < 16) is the last problem to be 
solved before a procedure can be devised for executing paths. 

One might sample by generating a string of exponential 
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random variables using the desired parameter lambda in order 

to get the time jumps < t . , t_, . . . , t > of the Poisson 

x £. n 

process which determines course change; then make a draw from 
each -from the x and y normal distributions at times t = tj, 
tj + t 2 »*-*> *n-l + *n *° ot) tain x ar| d y coordinates for the 
path legs. While the mean of the y coordinate draw is always 
zero, there is nothing in this method to prevent one y 
coordinate draw from falling on one side of the x axis and 
the next draw on the other side. As the cumulative time 
approachs T/2 the variance of the draw becomes relative 
large and it is quite likely that two consecutive draws which 
fall on opposite sides of the mean will be very far apart 
relative to the correspondi ng time step. A similar argument 
applies to the x coordinate, and because position is 
determined by both coordinates, if two consecutive draws are 
on the opposite sides of their means for both the x and y 
coordinates, the distance between two consecutive positions 
could be very far apart. The result is that a ridiculously 
high velocity is required in order for the target to travel 
from one sample point to its successor in the given time 
step. What is needed is some way to guarantee that a 
successive position is tied to the one before it. Future 
position must be conditioned on present position, and that is 
what is missing in this method. One possible remedy might be 
to reflect each draw across its mean if needed. Thus, the 
variance of the draw is preserved while the distance between 
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successive positions is shortened to a manageable length. 
However , this type of optional reflection does not capture 
the true Markovian property of each draw* but it does provide 
insight into how each successive draw should be done. 

Suppose that the first draw is carried out without 
reflection as described above. The first leg will connect the 
origin with a point selected by a draw from the two normal 
distributions < 16 ) substituting t = tj (the first time 
selected from the exponential distribution). Invoking the 
Markovian property, all information prior to this latest 
point is superfluous. Hence, there exists a completely new 
controlled time of arrival problem, wherein the latest point 
becomes the new origin and the line between that point and 
the destination becomes the new baseline. The total run time 
for the new problem is T - t ^ . This iterative process 
continues until the original total time <T) has expired. On 
the final leg, the cumulative time for the two coordinate 
draws is t = tj + t£ + ...+ t = T, so that the variance in 
(16) goes to zero and the target visits its destination on 
time. This procedure is a sensible one that seems to capture 
the properties of Brownian motion while at the same time 
producing discrete linear legs of target travel. Thus, it is 
the procedure adopted to produce the second type of path 
developed in this thesis, and the steps of the recipe are: 

1 . Select baseline length (distance between origin and 
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destination = L) , total nun time <T) , turning rate (\J> , 
maximum and minimum target speed, and physical 

parameter -for variance <d^> . Simulation frame of 

reference is a Cartesian coordinate system in 2-space, 
with origin at (0,0) and destination at <L,0). 

2. Generate a string of random variables from an 

exponential distribution with parameter lambda, and 
truncate the nth value in the string so that the sum of 

the random variables t, + t„ +...+ t = T. 

12 n 

3. Make a draw from the normal distributions in < 16) , 

using t = cumulative run time. For the first draw, 
t = tjj for the second draw, t = t^ + t 2 ; and so on. 

4. Measure distance between new point and previous 

point If the distance/t^ is greater than 

the maximum target sp>eed, truncate the leg as 

illustrated in Figure 3 if the distance/t. is less than 
the minimum target speed, extend the leg. Either 

truncation or extention results in new <X.,Y.). 

5. Reframe problem as a completely new one, using (X.,Y.) 

as the origin, <L,0> as the destination, and T = time 
remaining. Execute steps 3,4 again. The last draw 
forces <X,Y) to be equal to (L,0) unless truncation or 

extention occurs, in which case problem ends on time 

with target short of destination, or target continues 
to destination, in which case total travel time exceeds 
T. One path is completed. 
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Figure 3: 

Leg Truncation and Extension 



Figure 3a: Truncation 




I leg i | < min. spd. x t ^ ; 
so | le g i (new)| = min. spd. x t^ 
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6. Generate several hundred such paths and compute the 
statistics of interest -for each one at specific time 
intervals (mean square radial distance between present 
position and baseline position of original problem; and 
absolute difference between target present course and 
course to destination from present position). Look at 
plots of paths to determine if they * 1 ook good", as 
discussed in Chapter II. 

After statistics and plots are generated for both the 
random tour and Brownian motion derived paths, the task is to 
compare the two approaches in order to determine their 
respective strengths and weaknesses. 
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V. RESULTS 



A. RANDOi TOUR WITH DRIFT 

One measure of effectiveness for judging paths requires 
that they “look good", as detailed in Chapter II. It makes 
sense to try this measure on a path procedure first before 
expending effort to evaluate a mass of statistics. If the 
candidate paths can be rejected on sight then time will not 
be wasted on the other two measures. 

Figures 4a-f are represen tat i ve paths generated by the 

random tour with drift method. The circle in these figures, 

as described earlier, is inscribed with the end of the left 

path as its center and has a radius = V x time remaining in 

w 

problem after the execution of both the left and right paths. 
Thus, the target can easily travel, from the end of the left 
leg to the end of the right leg in the allotted time at the 
stated water speed. The ratio L' cJ /V^ listed on each of the 
figures measures just how constrained each path is. The 
ratio ranges from zero for the unconstrained case (random 
tour without drift) to positive infinity for the totally 
constrained case (straight line between origin and 
destination). Figure 4a shows a path for which ~ 2.69; 

clearly this path does not “look good* and is not acceptable. 
It is almost a straight line between origin and destination, 
and future target position is easily inferred from only a few 
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Figure 4c: Random Tour With Drift 
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Figure 4e: Random Tour With Drift 
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Figure ‘if: Random Tour With Drift 
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fixes. Of course, any method for devising a path under such 
constraints is bound to fail. A target cannot randomize its 
motion if it must travel directly to its destination. The 
rejection of the path in Figure 4a is not a rejection of the 
idea of random tour with drift, but only of the constraints. 
Even though 2.6? is in the low end of the range of V > any 
ratio greater than about 1.58 produces paths that are 
unaccep tabl e. In Figure 4b, where the ratio is 1.13, the 
path is still highly constrained but looks much better. The 
drift constraint really begins to loosen in Figure 4c, and 
in Figures 4d-f the ratio is 0.26. Here, the randomizing 
qualities of the random tour are evident. The paths all 
"look good", and whether they are executable by a target 
depends only upon what the time and distance units are. For 
instance, if the time unit is one half hour, Figure 4f has 
the target turning every 7.5 minutes on the average. This 
turn rate is not realistic for a large target. However, the 
path in Figure 4d has a mean turn rate of two per hour and is 
certainly executable. The main concern at this point has 
been put to rest; the paths in Figures 4c-f for which the 
time constraints are reasonable, look good enough to warrant 
further evaluation. 

The next measure of effectiveness against which to 
evaluate the random tour with drift is the distribution of 
the magnitude of difference between present course and course 
from present position to destination, measured at specified 
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times during path generation. This measure is discussed in 
Chapter II and is illustrated in Figure la. Recall now the 
procedure for generating each leg of the tour. The direction 
of the randomizing velocity (V ) is selected by making a draw 

P 

from a uniform distribution, and its magnitude is 
determined by multiplying the randomizing speed (determined 
from desired target speed through the water and drift, as in 
Figure 2) by the time length of the leg selected from an 
exponential (X) distribution. While the distribution of the 
direction of the randomizing velocity may be uniform, the 
distribution of the direction of the velocity through the 
water is not. The solid spokes of the left circle in Figure 
10a represent randomizing velocity and are equally spaced 
angularly. The dashed spoken of the right circle in Figure 
5a represent the water velocity after the appropriate drift 
velocity is added, and are not equally spaced angularly. The 
angle (6^) that any water velocity vector (V^) forms with the 
x axis is a function of the angle (0^) of the randomizing 
velocity, the drift speed, and the randomizing speed. In 
Figure 5b; 



r = V / X d - V . / X, 

r d ’ 

where the randomizing velocity required to cause the desired 
mean water velocity, as derived in (?) and (18), is, 

V = ( V 2 - u 2 )* 

r w d 
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because the mean angle between randomizing velocity and 
course from present position to destination is 90 degrees. 
The x and y components of the randomizing velocity are; 

x 0 * r cose 0 y 0 = r sin$ 0 

and the x and y components of the water velocity are; 



X 1 - *0 * d 



y l = V 



Now , 



tanSj = yj/Xj 



tanOj = y 0 / <x 0 + d> 



<<<V /X) sine.) / [<<V /X) cose a > + d3 5 
r 0 r 0 



= si nO 0 / Ccos0 0 + <V d /V p ) ] . 



And final 1 y , 



d, = arc tan <sin0- / Ccos0 Q + (U / V > 3 > < 17> 

1 0 0 d r 



Notice that the redi str i bu t i on of courses is not a function 
of lambda, but is strictly dependent on the ratio 

Even having the distribution of 0 0 and Knowing 0^ as a 
function of 9 0 , obtaining the density of 9^ by mathematical 
analysis is very difficult, but can be circumvented somewhat 
satisfactorily by using Monte Carlo simulation. It is a 
simple matter to draw several hundred random numbers from a 
uniform <0,2lO distribution and then tranform them with (17). 
This procedure was executed and the results appear in Table 1 
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Figure 5: 

Effect of Drift on Course Distribution 



Figure 5a: 

Redistribution of Courses After 
the Addition of Drift 




Figure '5b: 

Representation of Course 
Angular Shift 




54 



compared with the statistics obtained during path generation. 
Some significant observations ares 

1. The mean magnitude of the difference between present 

course and course to destination obtained from path 

statistics, and the associated standard deviation, 

agree very closely with the statistics obtained for 6^ 

by the simple simulation described above. This result 

was anticipated because the expected position of all 

experimental paths lies on the baseline even though 

very few actual positions fall there. However, and 

very importantly, Figures 6a, b, c and 7a indicate 

that there is a strict maximum deviation angle for any 

ratio. This is true only if the drift vector 

points directly from present position to destination 

all the time. If position is always on the baseline 

(where the mean position for the entire process lies) 

# 

then a maximum deviation will exist and will be arcsin 
, a fact which Figure 7a makes convincing. 
However, during travel the target deviates above and 
below the baseline, and though there is a maximum angle 

between the drift vector and water vector, to this 

angle must be added the depression or elevation angle 
between the horizontal and the destination. For 

example, if the target's x-axis distance from 

destination is 10 units, as is also the y— axis 
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then 



distance, and the target is above the baseiine, 

the angle between the horizontal line going through the 

target and the line going between target and 

destination is 45 degrees. Now, the next leg out of 

that position might -form the maximum angle 

Carcsin (V /V .) 3 between it and the drift direction, 

r a 

Then, the pointing deviation which we are interested in 
is the sum of those two angles = 45 degree's + 

arcsin . Nonetheless, the mean pointing 

deviation angles and their standard deviations are 
remarkably similar for the simulated paths and the 
simulation of f<$j) , where one set of statistics i <> for 
the true distribution of courses and the other is for 
the’ conditional distribution given present position 
equal to mean target position. 

2. No significant difference is noted between the mean 

differences for left unclipped, left clipped, and right 
paths. This is evidence that not only are left and 
right matching paths symmetrical in distribution, but 
more notably that clipping the left side by 

constraining it to match a right side path has no 
effect on the frequency with which the target points at 
i ts dest i nat i on . 

3. Lambda, as expected, has no effect on course 

distribution; only the ratio does. 
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Figures 6a-c show the distributions of simulated course 
deviations -for high, medium and low ratios. Notice that -for 
a high ratio, greater than one in this context and signifying 
tight constraints, the distribution loads on the higher 
values, while -for ratios less than one, signi-fying loose 
contraints, on the lower values. The break point occurs when 
the ratio is equal to one. While this loading, -first on one 
end o-f the distribution range and then on the other, seems 
odd at -first blush, there is a reasonable explanation -for it 
that Figure 7 helps to illustrate. For high ratios the 
tran-f ormat i on o-f to 0. causes a high proportion o-f the arc 
length o-f the circle described by the maximum to be 
subtended by the higher values o-f Figure 7a illustrates 
this occurence and also provides a graphical representation 
o-f why maximum values o-f course deviation are low for high 
ratios. As the ratio gets larger, the drift vector in the 
positive x direction accounts for an increasing amount of the 
water velocity, causing the circle described by the maximum 
magnitude of the randomizing velocity to grow smaller. 
Notice also that when K>. = V (Figure 7b) the maximum course 
deviation is 90 degrees and all feasible values are equally 
distributed because they subtend equal amounts of arc length. 
When the ratio is any greater than one, it is immediately 
feasible to have deviations as great as 180 degrees, the 
maximum possible. However, the higher deviation values 
subtend less arc length of the circle described by maximum V 
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FREQUENCY 



Figure 6a: 

Simulated Distribution of Course Deviations 
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Figure 6b: 



Simulated Distribution of Course Deviations 
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Figure 7 : 

Redistribution of Courses After 
the Addition of Drift 




Figure 7a: 



V 



V 



r 



Figure 7b: 



V 



V 



r 



F igure 



7c: 

V 

r 



62 



than do the smaller values, as illustrated in Figure 7c. 

Recall though that these distributions are not actual course 

deviation distributions determined -from generating actual 

paths, but are conditional course distributions given that 

present posi ton equals baseline position (which is expected 

position of target over all possible paths). However, and 

most importantly, the mean and standard deviation of the 

actual path course deviation distribution can be . determined 

using this simple inexpensive simulation once V. and V are 

d r 

specified. Thus, the target knows how much it tends to point 

at its destination or away from its origin. 

As detailed in Chapter II, the third measure of 

effectiveness is mean square radial distance between present 

position and baseline position. If the density of the course 

distribution after the application of drift [f<6j>3 had been 

2 

obtained, then EtR^ 3 could be easily calculated for any time 

making the appropriate substitutions into < 6) and <5). 

Again, however, the derivation of the density of is very 

difficult and not necessary because the process can be viewed 

as a random tour carried out at the randomizing velocity, 

merely shifted right or left by drift velocity prorated for 

cumulative run time. At any time, the "origin" of this 

process, viewed as a driftless random tour with velocity , 

2 

is the baseline position. This way, EIR^ 3 is calculated by 
substituting and t directly into (4). However, the 

resulting figures can only be compared validly to the 
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experimental statistics obtained for the unclipped left side; 
the requirement -for left and right sides to match adds a 
further constraint that is not captured by a process that is 
merely an unconstrained random tour appropr i atel y shifted for 
drift. Thus, the statistics for the matching paths must 
first be compared to each other in order to test for symmetry 
and then compared to the experimental statistics for the 
unclipped left side, which themselves have been compared to 
the figures obtained from (4). All these comparisons were 
made with the following significant results: 

1. Experimental statistics for both the left and right 

matching paths support the assertion that both sides 
are symmetrical. 

2. Experimental statistics for ihe unclipped left path 

mean square radial distance agree very closely with the 
analytical figures obtained tty substituting into 

<4> . Experimental figures tend to be one to three 

percent higher, but any difference that small is 
acceptable as sampling error. Mean square radial 
distance grows with time, and for a given time t is 
less for greater lambda, as expected. Table 2 
summarizes the statistics for the various time checks. 
'Cleft' denotes clipped left paths, 'left' denotes 
unclipped left paths, 'actual' means obtained from path 
simulation, and 'expect' means obtained from equation 
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< 4) by substitution of appropriate . The holes in 
the table occur because, at lambda equai to one or two, 
checks were not made -for all times listed. 
Nonetheless, more than enough statistics were obtained 
to provide valid results. 

3. The ratios o-f experimental c 1 i pped/ex per j. mental 

unclipped mean square radial distance were calculated 
■for the le-ft paths and are summarized in Table 3. 

Linear least squares regression o-f the ratios against 
time indicates that mean square radial distance 
reduction caused by clipping (requiring a match with 
the right side path) is independent of and lambda 
but dependent on time. The mean reduction is about twc 

t 

per cent per time unit with a standard deviation of 

approximately 1.2 per cent over the range for which the- 
regression was done. On this range a linear fit is 
quite good, but notice that after fifty time units the- 
linear reduction results in a mean square radial 
distance of zero, which is ridiculous because it is out 
of range for the regression. Clearly, the reduction 
caused by requiring the right and left paths to match 
must be calculated over the appropriate range. Though 
two per cent does not sound like much of a loss, it 

becomes quite significant after awhile, as will be 

demonstrated when the random tour with drift compared 
with Brownian motion. Nonetheless, it is clear that 
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Mean Square Radiai Distance from Baseline Position 
at Time t for Random Tour with Drift 



fs* m 

re m 
P- 

m rj 



O V 

co co 



w cv c mm cT'i/*' 

\ ~ r~ sr* — S 3 m on 



— • c: c, — — • 

m ’j 

o m a*! m co v -j 

♦ Q. rr f\ p\> sj CD m rg 

CO X • • • • • • • • 

rv c; vr v"* o c% cv cc 

i-* mm * r>. r**. 

^ cm rv, 

•3 liv ct'i vr m so -o 

-> -J-'i co v/i rj 

r*» '■«: c. co m rs v* v 

r» vn *r <n pj cs pj vf> 



c_ r» 
— O 



— * <— < mm mm 

r c co co ut it 

vP lT. — m — * 



VO VO 
S 3 CO 



u u 

*-4 U-. 

1 * O 



■ JU U xJ 

O 3 0 



-Q 

£ 



mm lO m 



' —« mm 



68 



Table 3: Ratio of Left Clipped/Left Undipped Mean Square 

Radial Distance at Time t for Random Tour with Drift 
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mean square radial distance -for any path, given dri-ft 
speed, randomizing speed, and lambda, is very 
pr edi c tab! e . 

B. BROWNIAN-DERIVED MOTION 

Again, the first measure of effectiveness addresses how 

representative paths produced by a method “look", and Figures 

8a-g are representative of discrete Brown i an-der i ved motion 

as executed by the second procedure delineated in Chapter IV. 

The paths in Figures 8a-d "look* acceptable; they show no 

pattern of regularity and seem to be executable, depending on 

the time and distance units selected. Figures 8a and 8b were 

generated from identical random numbers, V ^/V^ ratio, 

X = 1.0, and maximum and minimum water speeds. The only 

2 

difference between them is the physical parameter <? which is 
5.0 for Figure 8a and 10.0 for Figure 8b. The paths look 
quite similar but do not look different enough for one to 
speculate about possible differences for the other two 
measures of effectiveness. It is worth noting that for 
d 2 — 5.0 there are more extentions and less clips than than 
for Of 2 =10.0. This result is to be expected because a 
higher variance in the normal draws in (16) should produce 
successive positions that are farther apart than those 
produced with a lesser variance. When mean square radial 
distance is discussed later, the amount of clipping and 
extending which takes place becomes significant, generally 
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Figure 8a: PLOT OF DISCRETE BROWNIfiN-DERIVED MOTION 
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Figure Bb: PLOT OF DISCRETE BROWN I HN- DER I VED MOTION 
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rigure 8c: PLOT Of DISCRETE BROWN I RN~DER I VED MOTION 
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Figure 8d : PLOT Of DISCRETE BROWN 1 ON-DERIVED MOTION 
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Figure 8e: PLOT OF DISCRETE BROWN I RN-DER I VED MOTION 
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Figure ef: PLOT OF DISCRETE BROWN I ON-DERIVED MOTION 
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Figure 8g: PLOT OP DISCRETE BROWN I fiN~DER I VED MOTION 
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causing tne difference -for that statistic -for two different 
^2 to be reduced. 

Figures 8b and 8c were produced with similar inputs 
except -for lambda, which is 1.0 and 3.0 respectivel y. The 
difference between these paths is readily apparent, but again 
it is impossible -from one example -for each set of inputs to 
draw any general conclusions about the effect o-f turning rate 
on other- two measures o-f effectiveness. Figures 3d-g were 
generated from identical inputs with the exception o-f V 
which is 0.12, 0.23, 0.5, and 1.98 respectively. As the 
ratio becomes higher the paths straighten out as expected. 
The path in Figure 8g does not "look* acceptable, which is no 
surprise given :he strict contraints. Random tour paths with 
similarly high ratios su-f-fered the same straightening. 

Notice also, that as the ratio becomes high the process is 
more constrained but less likely to obey the constraints. 
The goal of the original controlled time o-f arrival problem 
was to assure that the target started -from an origin and 
visited the destination at the end o-f the stated time period. 
It is clear in all four of Figures 8d-g that the constraints 
are violated, and furthermore, that as the 'v> d /V r ratio 
increases the constraints are violated more. For any path 
that does not end at the destination there remains three 
options. First, it can be rejected outright. Secondly, it 
is possible that the requirement to reach the destination 
exactly is not as important as ending the path on time, in 
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which case, the path is acceptable if it is not too straight. 
Lastly, if the time constraint is not as important as 
ensuring that the target visits the exact destination, then 
more time can be given to the target with which to travel 
directly to destination if the distance is not great. The 
problem of control will be the subject of further comment 
when mean square radial distance from baseline position is 
examined. It suffices at this point to recognize that the 
problem exists and will probably affect the other two 
measures of effectiveness. 

The second measure of effectiveness quantifies how much a 
target points directly away from i:s origin and to its 
destination. For the random tour i * was found that the 
amount of pointing, as measured by the magnitude of the angle’ 
between the target's present course and the course to the 
destination from present position, was a function only of the 
ratio VV and not a function of time or turning rate 
lambda. Those results provide a starting point for the 
investigation of pointing for Browni an-der ived motion. 

Histograms of the deviation angle appear for high, medium, 
and low ratios for Br own i an-der ived motion in Figures 

9a-c , and Table 4 lists the statistics for three < 5 ^ values 
and six ratios. The table also lists values of pointing 
statistics for the random tour with drift obtained during 
path generation. Significant results are: 
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1. The amount of pointing for Br own i an -derived motion 

agrees very closely with that for the random tour with 
drift for low ratios, as exhibited in Table -4. 

This result is expected because, regardless of tne 
method chosen for generating paths, a drift velocity 
will bias course headings in the direction of the 
destination. Brown i an -der ived motion with the same 

V ratio as a random tour should exhibit similar 
pointing behaviour. 

2. Table 4 also shews that as the dr i f t/randomi z i ng ratio 

increases the pointing statistics for the two different 
processes begin to diverge. While the random tour 
statistics continue to agree closely with those 
generated by the simulation of (17), the Brownian- 
derived motion statistics decrease to a mean of 

approx imatel y 5 6 degrees and standard deviation of 
about 20 degrees. Here, the effects of loss of control 
at high dr i f t/r andomi z i ng ratios appear. The pointing 
statistics are more favorable, but at the cost of the 
target breaking constraints. Again, imagine the stick 
with a piece of string attached to each end as 
described in Chapter II. If one end of the string is 
loosened and required only to be near the end of the 
stick, that is equivalent to adding more string and 
thereby easing the constraints. When the mean square 
radial distance for Brownian-derived motion is 
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examined, the lack of control at high drift/randomizing 
ratios will become evident. 

3. The pointing statistics show that for Brown i an-der i ved 
motion lower values of <f^ produce slightly more 
pointing. This phenomenon is most pronounced at high 
V d /'v‘ r ratios, where the process is not under control 
anyway. However, it demonstrates that even at the same 
ratio, a smaller variance for the position sampling 
causes the process to be more constrained. The slight 
difference in statistics would most likely be greater 
if low variance legs were not extended more and clipped 
less than high variance ones. 

4. The histograms in Figures ?a-c shew that the standard 

deviations of the distibutions tend to increase and 
decrease with the means. The low mean distributions 
are peaki«-r than the high mean ones, indicating 

processes that are more constrained; the opposite is 
true of less constrained paths. 

The third and final measure of effectiveness for Brownian- 
derived motion is the mean square radial distance between present 
position and baseline position, measured periodically throughout 
path generation so that a graph or set of ordered pairs 
<t, * s produced for comparison with other path generation 

methods. Significant results for this measure are: 



82 
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FRFQUENC 1£S SAMPLE SIZE * 7250 
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1. Figures 19a-o illustrate mean square radial distance 

■for four V ^/K> ratios and X = 1.9; each figure has 
curves for three different values of the physical 
parameter d Immediately it is clear that the mean 
square radial distance decreases as the 

dr i f t/r anoomi 2 i ng ratio increases. This result follows 
directly from the greater constraint placed on the 
target at higher ratios, and is similar to the results 
for the random tour with drift. 

2. The problem of control for Brcwnian-derived motion is 
apparent in these figures. Even in Figures 10a and 
10b, at the lowest two dr i f t/randomi z i ng ratios, the 
mean square radial distance does not go down to zero, 
indicating that on the average the target does not 
visit the destination. Notice that as the ratio becomes 
progressi vel y higher in Figures 10c and 10d that the 
mean square radial distance tends to become even 
greater at the end of the time period. Clearly, as the 
problem becomes more constrained, control is lost. 
Another characteristic of this path generating method 
with small lambda is the hook in the curve at the end 
of the time period, indicating an undesireable increase 
in mean square radial distance during the final time 
unit of travel. This increase is the result of the 
truncation of the final leg when time runs out, though 
the random numbers for the final leg were generated as 
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Figure 10b: MEAN SQUARE RADIAL DISTANCE 
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Figure lod: MEAN SQUARE RADIAL DISTANCE 
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TIME 



if no truncation takes place. The simple solution is 
for the target to travel directly coward the 
destination on the last leg in order to get as close as 
possible. A target would probably do that anyway since 
it does not really need to be told what to do on the 
final leg. The algorithm for path generation was 
purposely not modified in order to demonstrate this 
peculiar phenomenon. However, it is important to 
realize that a true "bridge", Brownian or otherwise, 
does not exist as long as the two endpoints are not 
connected by the target path. 

3. At lower drift/randomizing ratios, the mean square 
radial distance is less for smaller as Figure 10a 

clearly illustrates. However, as the ratio increase's, 
the mean square radial distance is approx imatel y the 
same for all three values of <?£• The curves in Figures 
10b-d also show that for different d 0 values the 
maximum mean square radial distance occurs at different 
times for a given dr i f t/randomi z i ng ratio. The mean 
square radial distances for different < 3 ^ become closer 
in value because of the clipping and entending process 
that occurs to keep the target within speed limits. As 
the problem becomes more constrained at higher drift/ 
randomizing ratios, the target tends to bump into the 
speed bounds more. Hence, the process is clipped and 
extended more and the natural differences among the < 3 ^ 
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which all 



values dominate less than the speed limits, 
paths must obey equally. Also, as the mean square 
radial distances -for various sigma square all get 
closer to the baseline for whatever reason, they are 
bounded on one side by the baseline (the minimum value) 
and get sandwiched together, as Figure 19d clearly 
demonstrates. 

4. Figures 10a-d show that for a given V c j/^ p there is a 

value of sigma square which produces a mean square 
radial distance curve that is symmetrical about the 
line x = 1 / 2 . For Figures 10a-c the values of sigma 
square which come closest to symmetry are 19, 5, and 1 

respec t i ve 1 y . Notice also, that in Figures 10a and 

10b, in which the dr i f t/randomi z i ng ratios are low and 
sigma square equal to ten and five are the values 
closest to producing symmetry, the curve for sigma 
square equal to one is very asymmetr i cal and flattens 
out for big t. This occurs because for these low 
dr i f t/randomi z i ng ratios, = 1 over-constrai ns the 

process; it forces the target back to the baseline too 
soon . 

5. Figures lla-d are different from Figures 10a-d only in 
their lambda value which is 3.0 instead of 1.9. It is 
immediately clear that an increase in lambda causes a 
marked decrease in mean square radial distance. A 
similar result was found for the random tour, for which 
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an increase in lambda caused the distribution of 
radial distance to start piling up around the expected 
position or the process. There are some other subtle 
di f f er ences also. At low ratios, all curves for 
the three sigma square values are skewed slightly less 
left than for higher lambda, and for higher 
dr i f t/randoini 2 i ng ratios they are skewed slightly more 
right. Notice also that the curves for the lowest 
dr i f t/r andomi 2 i ng ratio in Figure 11a are closer 
together than they were for the same ratio at 
X = 1 . <i in Figure 10a. This, again, is the result of 
all curves being bound on one side by the baseline, and 
therefore forced together more at the lower mean square 
radial distances* which higher lambda produces. 
Final tv, the "hook" on the end of the curves is gone, 
most likely because the mean number of turns during the 
final lime unit is three instead of one. Consequently, 
the target turns a few more times during the last time 
unit of travel toward the destination and thus does not 
rely only on one leg to hit or miss. Again, however, 
it is likely that a target would travel directly to 
destination at this point anyway, unless it desired to 
adhere to its course change policy strictly. 
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Figure lid: MEAN SQUARE RADIAL DISTANCE 
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CONCLUSION S 



Chapter I defined a controlled time of arrival scenario 
in which a time constraint is placed on a target that is 
required to travel between two endpoints. It is desireable 
•for the target to “randomize 0 its motion during the transit 
in order to provide the enemy with 3S little in-formation as 
possible about target origin and destination and to make 
target detection and redetection di-f-ficult. Chapter II 
discussed desireable qualities -for such target travel and 
delineated three measures o-f e-f -f ec t i veness against which to 
measure any procedure -for producing target paths, and 
Chapter III examined previous investigations into the problem 
which have provided direction -for the two approaches adopted 
in this thesis. Chapter IV described the random tour with 
drift and discrete Brown i an -dervi ed motion in detail, and 
while the results for each method were presented in 
Chapter V, the performance of a path producing procedure 
against the measures of effectiveness is not important in 
itself; but rather, the measures provide a way to compare two 
or more procedures, one of which can be judged best for a 
particular situation. 

The amount a target points away from orign or to 
destination is very nearly the same for both the random tour 
with drift and discrete Browni an-deri ved motion. Tables 1 
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and 4, which list the mean values and standard dev*, a 1 1 on s of 
the course deviation distributions, show nearly identical 
values for both procedures at low ratios. Differences 
begin to show only as the ratio gets bigger. But though the 
statistics become more favorable for Brown i an -der ived motion, 
it is probabl y because that process is not under complete 
control of the constraints when the dr i f t/r andomi z i ng ratio 
is high. If it were completely controlled or if the random 
tour were also allowed to violate the final visitation 
constraint, then the statistics would probably be very nearly 
identical over the whole range of U ratios. 

The failure of Br own i an -der ived motion to control target 
travel by guaranteeing visitation on time at the destination 
is a potential weakness of the process; it literally does not 
do exactly what it is supposed to do. However, the weakness 
is not important just as long as it is not vital for the 
target to get exac 1 1 v to destination exactly on time. If it 
is vital, then clearly the choice is in favor of the random 
tour, and none of the measures of effectiveness are relevant, 
except for the requirement that the path “look good" and be 
executable as described in Chapter II. 

If the target does not need to have a perfectly 
controlled time of arrival, and can either fall short of its 
destination at the appointed time or else take extra time to 
get there, then the mean square radial distance becomes the 
important measure of effectiveness; pointing to origin and 
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destination is similar and representative paths of both 
procedures " 1 ooK good" and are executable. Figures 12a-d 
illustrate the mean square radial distance between present 
and baseline positions -for two sigma square values of 
Brown i an -der ived motion and for the random tour with drift 
with and without the requirement for the left and right paths 
to match. Nhile the mean square radial distances for random 
tour paths which do not match have no practical significance 
since such paths meet none of the constraints, they do show 
hew much radial distance is lost by requiring that the left 
and right paths meet. Recall that the loss is about two per 
cent per time unit, and while that sounds low, it is becomes 
very significant as time progresses. The figure of two 
percent was obtained, -as previously described, by linear 
regression and is valid only over the range of time for which 
the regression was done. This limitation is illustrated in 
Figures 12a and 12b. Notice that the curve for matching 
random tour paths begins to sag in the middle, when actually 
the mean square radial distance should be increasing, however 
slightly. The sagging is a direct result of the linear 
regression operating at the edge of its valid range. The 
curve should be rather flatly rounded in the middle instead 
of sagging. However , the curves presented are accurate over 
their range, with the exception of the slight sagging as 
described. In Figures 12c and 12d, where the time period is 
shorter, the curves do not sag at all. Most importantly, all 
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Figure 12a: MERN SQUARE RRDIRL DISTANCE 
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Figure 12b: MEAN SQUARE RADIAL DISTANCE 
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Figure 12c: MEAN SQUARE RADIAL DISTANCE 
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four figures show how much mean square radial distance is 
lost by requiring the paths to match. Under the conditions 
prevailing in Figure 12a, the unconstrained random tour with 
drift exhibits a greater mean square radial distance than 
Brownian-derived motion does *o r either sigma square value. 
However, the mean square radial distances for the matching 
random tour paths are much less than those for Brownian- 
dervived motion in all figures. Notice also that while an 
increase in the V ratio or lambda each cause a reduction 
in mean square radial distance for both procedures, that 
increasing lambda affects the random tour process much more 
adversely than Brown i an-der i ved motion. This is a point in 
favor of the Brown i an-der ived process. 

One might argue that it is natural for Brown i an-der i ved 
motion to exhibit higher mean square radial distances, if 
only because the process is not under control. This argument 
is not compelling and one need only to look at the figures to 
see the great disparity. It is not likely that the failure 
of the mean square radial distance for Brownian motion to go 
all the way to zero at the destination is the reason that it 
is more than double that of the random tour in the midrange 
for three of the four cases. However, one might argue 
further that the constraint which requires the left and right 
random tour paths to meet up causes such severe degradation, 
and that Brownian motion might suffer similarly if it could 
be made to meet the constraints exactly. Nonetheless, the 
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paths perform di f f er en 1 1 y as they are. and clearly the mean 
square radial distances are much more favorable for Brownian- 
derived motion if the strict breaking of the final visitation 
constraint can be tolerated, and remedied. 

Recall that the first measure of effectiveness which is 
applied to a path producing procedure checks representative 
paths to see whether they "look good* and are executable by a 
candidate target. Failure against this measure au tcma t i cal 1 y 
disqualifies a path from further consideration. While both 
the random tour with drift and Brownian-derived motion passed 
this vital first test, there is a subtlety in the way both 
procedures are executed that makes Brown i an-der i ved motion 
more desireable. In order to obtain matching left an d r: ght 
paths, the random tour must be executed in its entirety 
before the target begins a journey. The target must "hen 
follow the instructions closely to make all the correct 
courses and turns. While this procedure is possible to carry 
out, it is quite exacting. On the other hand, Browmian- 
derived motion can be executed one leg at a time because 
after each leg the controlled time of arrival problem is 
reframed as a totally new one using present position as the 
new origin. Hence, perfect navigation is not as critical as 
it is for the random tour. In a sense, starting the process 
over after each leg always gives the target another chance, 
just so long as a flagrant violation does not occur which 
causes the target to be faced with an impossible transit at 
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the beginning of some intermediate leg. This -feature of 
Brown i an-der ived motion should make it more saleable to 
target captains. 

Both methods -for producing paths should be available -for 
a target to choose. Clearly, i-f the -final visitation 
constraint absolutely must be met on time, then the random 
tour with dri-ft method provides the only guarantee. But, i-f 
the target needs only to ensure that it arrives in the 
immediate vicinity o-f the destination or may arrive at the 
actual destination slightly ear 1 y or late, then discrete 
Brown i an-der ived motion per-forms more -favorably against the 
selected measures o-f e-f -f ec t i veness and is easier -for a target 
to execute. 

The two very di-f-ferent methods presented here -for 

generating paths represent only two among many, and the 
variations on these two procedures alones are infinite. For 
instance, the random tour with drift could be modified so 
that it was executed from one end only, instead of from both 
ends as done in this thesis. After each leg of travel the 
problem could be reframed as a totally new one, in much the 
same manner that Brown i an-der i ved motion was restarted after 
each leg here. It is also quite possible, and desireable, to 
devise a way to force Brown i an-der ived motion to arrive at 
the destination exactly on time, satisfying all the 

constraints strictly. Thus, the two methods presented here 
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are not only quite workable, but also provide 
■for -further investigation. 



suggest ions 
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